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I SUMMARY 

^ Knowing the stress and strain distributions in the left ventricular wall 
jf the heart is a prerequisite for the determination of the muscle elasticity 
'nd contractility in the process of assessing the functional status of the 
leart. NASTRAN has been applied for the calculation of these stresses and 
[trains and to help in verifying the results obtained by the computer pro- 
jram FEAMPS which was specifically designed for the plane-strain finite-element 
na lysis of the left ventricular cross sections. Adopted for the analysis are 
he true shape and dimensions of the cross sections reconstructed from multi- 
'lanar x-ray views of a left ventricle which was surgically isolated from a 
log’s heart but metabol ica I ly supported to sustain its beating. A preprocessor 
[as been prepared to accommodate both FEAMPS and NASTRAN, and it has also 
facilitated the application of both the triangular element and Isoparametric 
|uadr i latera I element versions of NASTRAN. The stresses in several crucial 
Regions of the left ventricular wall calculated by these two Independently 
level oped computer programs are found to be in very good agreement. Such 
■onfirmation of the results is essential in the development of a method which 
[issesses the heart performance. 


'’This work was supported in part by research grants HL4664, HL3552, and RR7 
from the National Institutes of Health, and by the Engineering Research Center, 
University of Nebraska-L i ncol n . 
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INTRODUCTION 


I 


Structurally speaking, human and animal hearts are very complex not only 
in geometry but also in material characteristics. The finite element method 
therefore presents itself as an effective tool for analyzing the stresses and 
strains in the hearts. In order to assess whether or not a human or animal ^ j 
heart is functioning normally, it is necessary that the working characteristic 
of its left ventricular wall muscle (myocardium) be adequately ascertained, fd 
it is the contraction of the left ventricular myocardium that pumps the 
oxygenated blood to circulate in the body. For the investigation of the left 
ventricle of the heart alone (not to be involved with the other chambers, righ 
ventricle, and left and right atria), experiments of surgically isolated but 
metabol leal ly supported, beating left ventricles of dogs are often prepared so 
that the dynamic changes of the shape and dimensions of the cross sections of 
fhe left ventricles during cardiac cycles can be recorded on videotape and by 
use of reconstruction techniques involving multiplanar x-ray projections or 
ultrasonic echoes of the cross sections (ref. I). 

The reconstruction technique using x-ray is schematically illustrated in 
Fig. I. For a certain cross section of interest, a number of x-ray projection 
are taken by rotating the cross section with a chosen angle increment about it 
normal axis. The shape and dimensions of the cross section can then be recons 
tructed by use of algebraic algorithms (ref. 2). With a life supporting syste 
(Fig. 2) for the isolated left ventricle of a dog experiment, any desired 
anatomic site selected at 0.7 mm intervals along the entire ap ica |-to-basa 1 ax 
of the ventricle can be reconstructed. Figure 3 shows a sample rectangular 
array of 35 ap ica l-to-basa I cross sections at a certain Instant of a cardiac 
cycle. Further details and discussion of the reconstruction and data collect! 
techniques are given in an Internal publication (ret. 3), which Is available 
upon request. 

During relaxation (diastolic) phase of a cardiac cycle, the myocardium 
extends as the blood fills into the left ventricle. Based on the values of th 
left ventricular pressure and volume concomitantly recorded during the dog 
experiment, the history of the external work being done to the left ventricle 
can be calculated. Meanwhile, the finite element method can be applied for th 
analysis of the reconstructed instantaneous shapes of the left ventricular cre 
sections to determine the stresses and strains in the wall muscle and conse- 
quently the change of the strain energies. Upon relating the external work ar 
the strain energies, the dependency of the elastic stiffness of the myocardiun 
on the chamber pressure can then be estimated. 

While the calculation of the external work is straightforward, it is for 
the calculation of the stresses and strains in the left ventricular cross 
sections that this paper reports the application of NASTRAN as a verifier for 
the plane-strain finite-element computer program FEAMPS which was specificallv 
developed for the study of the left ventricles. 
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LEFT VENTRICLE MODELS AND FEAMPS 


' Sine© the thin—walled shell model was found to be adequate only for calcu~ 
ating the mean stresses across the wall thickness of the left ventricle, a 
umber of thick-walled models have been proposed for obtaining the variation 
!f the wall stresses (ref. 4). Most of the thick-walled models approximated 
'he geometry of the left ventricle as a shell of revolution, either as a circu- 
ar hollow cylinder, hollow sphere, or hollow prolate spheroid, all having uni- 
form wall thickness. When biplane silhouettes of the left ventricle could be 
;aken by application of the roentgen-, cine- and videodensitometry techniques 
ref. 5), it became possible to incorporate some of the dimensions, such as the 
ariations of wall thickness and chamber radii, of the left ventricle into the 
hick-walled models to determine the wall stresses by axisymmetric finite- 
llement analysis. Now the advent of cross-sectional reconstruction from 
lultiplanar x-ray views enables the thick-walled models of the left ventricle 
[o be further improved by adopting the true shape and dimensions of the left 
fentricular cross sections. 

' However, the 35 apical -to-basa I cross sections shown in Fig. 3 were ob- 
ained by a reconstruction system using only one x-ray source and one detector, 
he system requires the position and geometry of the left ventricle to be kept 
onstant throughout successive cardiac cycles. This "stationari ty" requirement 
,s a serious disadvantage to be overcome by use of a spatial reconstruction 
lystem which will have multiple x-ray sources and multiple detectors. At 
Present, the available data include (I) the dynamic changes of the cross- 
lectional shapes and dimensions at a few sites along the apical -to-basa I axis 
k the left ventricle but are complete for several cardiac cycles, and (2) the 
|hapes and dimensions for a I Lapical-to-basal cross sections of the left ven- 
tricle but only at a few crucial instants of cardiac cycles. So these data 
fUffice to investigate the cross-sectional behavior but not the ultimate three- 
Hmensional analysis of the left ventricle. 

I FEAMPS is an abbreviation for Finite-Element Analysis of Myocardium by 
^lane-Strain Theory. The left ventricle is approximated as a hollow thick- 
/alled cylinder with uniform cross section having the shape and dimensions re- 
constructed from multiplanar x-ray views. When the cylinder is subjected to 
jniform chamber pressure, it is assumed that all cross sections will deform 
identically. By looking at the silhouette shown in Fig. 4 constructed from 
stacking all cross sections in Fig. 3, it appears that plane-strain analysis 
should give a reasonably good prediction of the cross-sectional- behavior in 
I'he vicinity of the base of the left ventricle. 

[ 

It has been well established that the left ventricular myocardium is 
fibrous. Several investigations of the left ventricular wall with fibrous mate- 
■-ial properties had been reported (refs. 6,7). However, only scattered data are 
available regarding the changes of the fiber directions during cardiac cycles; 
fhe results of these investigations remain to be verified by more thorough ex- 
aeriments. FEAMPS follows the approach of assuming the left ventricle as a 
lomogeneous, isotropic and elastic medium and proceeds to evaluate the effec- 
Nve elastic stiffness for such an equivalent model (refs. 8,9). Since the 
affective elastic stiffness is to be inversely determined on the basis of the 
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work-energy principle described in the Introduction, its value is assumed to 
be equal to I In the Incremental -loading analysis of FEAMPS. Regarding the 
other material constant, Poisson's ratio, of the left ventricular myocardium, 
the Issue whether it should be assumed equal to 0.415 (ref. 10) or nearly 
Incompressible equal to 0.49 (ref. ID Is still not yet settled. So the re- 
sults for both values will be discussed. 


FORMATION AND ANALYSIS OF ELEMENTS 


Figures 5(a) and 5(b) show a typical change of the shape and dimensions oi 
a cross section near the base of the left ventricle reconstructed from multi - 
planar x-ray views during the early and ending stages of the diastole of a car- 
diac cycle, respectively. The diastolic phase Is the resting period of the hed 
during which the relaxed left ventricular chamber is filling with blood and 
the wall muscles are passive and extended. The two-dimensional In-plane trl- | 
angular elements of FEAMPS are generated by adopting 9 nodes across the wall 
thickness and 30 nodes around the perimeter of a reconstructed cross section. , 
Additional elements may be necessary for the papillary muscle at certain levels 
within the ventricle. Figure 5(c) Illustrates that the end-dlastol Ic cross 
section shown in Fig. 5(b) is partitioned into 508 triangular elements using 
290 nodes, of which 20 are for the formation of the 28 triangular elements for 
the peninsulalike papillary muscle. 

The only input data required for FEAMPS are the coordinates of 30 pairs o1 
points on the inner (endocardial) and outer (eplcardlal) boundaries of the 
cross section, plus additional pairs if necessary for the papillary muscle. A 
preprocessor which consists of several subroutines has been prepared to carry 
out the sequential numbering and calculation of the coordinates of all nodes, 
to compose the interconnected triangular elements by defining their respective 
constituent nodes, and to generate the nodal external loads. Plot subroutines 
are also made aval lable for drawing of only the borders of the cross section 
or the partitioning pattern of the cross section either with or without the 
labelling of the node and element numbers. 

The plane-strain finite-element analysis of FEAMPS follows closely the prc 
cedures delineated in Zienkiewicz’s book (ref. 12). The displacements at any 
point within the triangular element are assumed as linear polynomials of its 
coordinates. This results in constant strains throughout the element. For 
saving computer storage spaces, the structural stiffness matrix is formulated 
rectangular form by taking advantage of symmetry and discarding a I I of the zer< 
elements that are on the outside of the diagonal band. In order to handle suci 
a rectangular coefficient matrix, the Gauss-Jordan elimination method has been 
accordingly modified for the solution of the nodal displacements. This approai 
has reduced the computer storage core requirement for the structural stiffness 
matrix from 580x580 to 580x22 floating words in the analysis of the end- 
diastolic cross section partitioned as shown in Fig. 5(c). In case that further 
reduction of core storage in the computer is necessary, the zero elements inside 
the diagonal band of the structural stiffness matrix may also be discarded and 
all nonzero elements compacted into a vector, and the Gauss-Seidel iteration 
method can be applied for the solution of the nodal displacements (ref. 13). 
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, Since the left ventricular cross sections are subjected to uniform internal 
hessure only on the inner (endocardial) boundary and there are no boundary con- 
fraints at all, the nodal displacements determined by machine computation may 
Dntain the superposition of a free body motion. To ascertain whether the some- 
Jmes very large free body motion will introduce significant errors in the sub- 
bquent calculation of the cross-sectional stresses, studies of the cross 
action with added fictitious boundary conditions have been conducted. More 
stalls on these studies are presented in the Discussion of Numerical Results. 

FEAI'IPS also has the built-in provisions for the calculation of the nodal 
tresses by averaging the stresses of the adjoining elements, the extrapolation 
f the stresses for improved accuracy at the boundary nodes, the calculation of 
he strain energies, and for various plottings of the stress distributions of 
he cross section. Figure 6 displays a sample plot of the nodal circumferential 
tress distribution (normalized with respect to the uniform internal pressure) 

,f the end-diastolic cross section with the exclusion of the almost uniformly 
pmpressed papillary muscle. Full details and computer program listing of FEAMPS 
re given In ref. 3. Nevertheless, the above synopsis should suffice for the 
[iscusslon of the present paper. 

1 

APPLICATION OF NASTRAN 


I 

Attempts have been made to apply NASTRAN Rigid Format I to the two- 
jimensional plane-strain static analysis of the reconstructed left ventricular 
iross sections. With the availability of the FEAMPS preprocessor, preparatory 
orks that can be expedited are the generation of the coordinates of the nodes, 
ihe division of the cross section into either triangular TRMEM or quadrilateral 
pMEMI elements, and the conversion of the uniform internal pressure on the 
Inner (endocardial) boundary of the cross section into equivalent concentrated 
oads at the endocardial nodes. 

I As has been mentioned in the preceding section, there are no boundary con- 
itraints at all in the plane-strain analysis of the left ventricular cross 
sections. Test runs are necessary for the determination of acceptable ficti- 
ious points of support which can be adopted as the constraints data for the 
5PCI cards and will not introduce additional external loading. For instance, 
in the analysis of the cross section shown in Fig. 5(c), the decision to con- 
strain node 6 in both the x and y directions and node 19 in the y direction 
/as based on the verification by the printout, via SPCFORCE request, that the 
^eactive forces at these two nodes were insignificantly small relative to the 
External loads at the endocardial nodes. 

Since the left ventricle is assumed as a homogeneous, isotropic and elas- 
Nc medium, the f^T2 cards specifying the material properties were prepared 
vith the values of the Young's modulus, E, equal to I and the Poisson's ratio, 
equal to 0.415 and 0.49, same as for the FEAMPS cases. 

The extreme stresses are known to occur along the inner boundary of the 
left ventricular cross section; the four regions of primary concern have been 
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particularly pointed out in Fig. 7 to facilitate the comparison study of the 
FEAMPS and NASTRAN results. These regions are the inner borders of the pos- i 
terior, septal, anterior and free walls in the vicinities of nodes 247, 253, 

258 and 264, respectively. Figure 8 gives a close-up view of these four 
regions with the node and element numbers indicated. i 


DISCUSSION OF NUMERICAL RESULTS 


Level 15.5 was used to generate all NASTRAN results reported herein. 

Table I shows that for v=0.4l5, the values of the principal stresses (normalizt 
with respect to the chamber pressure) in the end-diastolic cross section deterj 
mined by FEAMPS and I^ASTRAN are almost identical when both used triangular elei 
ments. The last two columns of Table I reveal that the principal stresses in 
quadrilateral element are very close to the average values of those of the two 
adjacent triangular elements which have been combined to form the quadrilateral 
element. Since the number of the elements is reduced to half for the QDMEMI j 
analysis but the computer time is about doubled, the TRMEM analysis was selects 
for further studies. 

For v=0.4I5 but for the early diastolic cross section shown in Fig. 5(a), 
Table 2 is presented to depict the effect of imposing fictitious boundary con- 
straints to the cross section on the principal stresses. It is clearly indica- 
ted that restricting the x-d i sp lacement u of node 6 and the y-d i sp I acement v 
of nodes 6 and 19 causes no significant change to the values of the stresses. 

It also helps to demonstrate that the free body motion of the cross section 
which may result by FEAMPS analysis will not alter the values of the stresses. 

The cases of treating the left ventricle as a nearly incompressible mate- 
rial with v=0.49 have also been investigated. Table 3 summarizes the norma I i zi 
principal stresses for both the early and end diastolic cross sections shown ii 
Figs. 5(a) and 5(b). The results exemplify the changes of the stresses during 

the diastolic phase of a cardiac cycle. As in Tables I and 2, these results 
obtained by FEAMPS and NASTRAN are in very good agreement. 

Since FEAMPS has a feature for calculating rhe nodal stresses, the normal- 
ized principal stress distributions through the wall thickness at the four re- 
gions of concern during early and end diastole have been obtained and are pre- 
sented in Fig. 9. It is noteworthy that some compressive circumferential 
stresses may develop at the innermost ( 0 % of wall thickness in Fig. 9) myo- 
cardial layer in the septal and free walls where the curvature are small and 
sometimes even turn convex inward. Because other axi symmetric models of the 
left ventricle assumed cross sections to be of annular shape, in no portion of 
the cross section could the wall be curved inward. As a consequence, no 
compressive circumferential stress in the left ventricular wall had ever been 
determined by these models. 
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CONCLUDING REMARKS 


I The evaluation of the cross-sectional elastic stiffness (ref. 14) of the 
left ventricle during diastole depends in large measure on the accuracy of the 
yalues of the stresses and strains in the cross section. So it Is indeed re- 
assuring that these values can be confirmed by FEAMPS and NASTRAN. 

This comparison study has also proved the applicability of NASTRAN for the 
plane-strain analysis of the left ventricular cross sections. Certainly, the 
other capabilities of NASTRAN can further help resolve some of the problems 
connected with the ongoing research of the three-dimensional analysis of the 
left ventricle, for which a special-purpose computer program is being developed. 
In the computer program, the true three-dimensional geometry of the left ven- 
tricle determined by the multiplanar x-ray views of its cross sections will be 
utilized and the analysis by treating the myocardium as a fibrous material using 
adjusted elastic stiffnesses and incremental loading will be incorporated. 

It may also be anticipated that NASTRAN will contribute to the ultimate dy- 
namic analysis of the left ventricle which will consider, among others, the ef- 
fects of the myocardial contraction, the dynamics of the heart valves, and the 
iinteractlon of the blood flow in the left ventricular chamber. 
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Table 1. Comparison of principal stresses for the end-diastolic cross 


section (E«p^«l and v«0.415). 
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Table 2. Comparison of principal stresses for the early diastolic cross i 


section (E*p^»l and v»0.415). 


I 

Element 

NASTW 

iN 

FEAMPS 1 

'l9*0 


— 

m\ 

Number 

*^1 


^1 

®2 

‘^i 

°2 

H 


BBS 

1.44 

0.48 

m 

1.44 

0.48 

5.01 

2.76 

1.44 

0.48 

Hi 

mmMl 

^KEEH 

1.96 

7.14 

0.30 

2.73 

1.96 

7.14 

0.30 

2.73 

1.97 

7.14 

0.30 

2.73 



5.51 

0.15 

0.06 

-1.38 

5.51 

0.15 

0.06 

-1.38 

5.51 

0.15 

0.06 

-1.38 



1.66 

2.63 

-0.31 

0.11 

1.66 

2.63 

-0.31 

0.11 

1.66 

2.63 

-0.31 

0.11 

1 


0.44 

-0.10 

-0.83 

-1.21 

0.44 

-0.10 

m 

0.44 

-0.10 

-0.83 

-1.21 

^KE^H 

-0.11 

-0.81 

m 

-0.11 

-0.81 

s 

-0.11 

-0.81 

BKIW 

WoEm 

— 

0.004 

-0.01 

m 

0.005 

-0.005 

m 

o.oor^ 

-0.01 

s 


s 

m 



I^^E 

s 



0.11 1 
2.07 

-1.59 

-0.01 

0.11 

2.08 

ESI 

Eb9 

0.11 

2.08 

-1.59 

-0.10 

448 

449 

1.51 

2.63 

-0.97 

-0.17 

1.51 

2.64 

-0.97 

-0.17 

1.51 

2.64 

-0.97 

-0.17 


3.13 

3.60 

kbm 

3.13 

3.60 

m 

3.13 

3.60 

-0.96 

-0.49 

H 

2.45 

2.54 

-0.56 

-0.50 

2.46 

2.55 

s 

2.46 

2.55 

-0.56 

-0.50 

F 

R 

E 

E 


4.08 

2.40 

m 

4.09 

2.40 

s 

4.09 

2.40 

mm 

456 

457 

1.80 

0.99 

-0.56 

-1.19 

m 

m 

1.80 

0.99 

-0.56 

-1.19 


1.71 

0.82 

m 

1.71 

0.82 

m 

1.71 

0.82 

-0.22 

-1.22 


0.72 

0.58 

-0.88 

-1.07 

m 

-0.88 

-1.07 

0.72 

0.58 

-0.88 

-1.07 
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Table 3. Comparison of principal stresses for the early and end diastolic 

® nearly incompressible material 

v=0.49 (E=p^=l). 


R 

E 

EARLY D 

lASTOLE 

END DIASTOLE 

G 

I 

n 

FEAf 

IPS 

NASTRAN 

FEAMPS 

NASTRAN 

u 

N 

®1 

®2 


‘^2 

®1 

^^2 

*^1 


P 

5.69 

2.81 




0.26 


0.19 

0 

3.87 

1.86 


M.M 


-3.59 


-3.45 

S 

1.83 

-0.08 

1.83 

-0.09 

1.63 

-1.40 

1.55 

-1.37 

T 

r* 

8.40 

4.65 

8.41 

4.65 

-3.25 

-7.33 

-3.14 

-7.06 

E 

p 

6.96 

0.29 

6.98 

0.29 

15.23 

10.21 

14.56 

9.74 

i\ 

I 

-1 .54 

-3.42 

-1.53 

-3.42 

-0.15 

-4.69 

-0.15 

-4.50 

0 

2.72 

0.68 

2.72 

0.67 

10.31 

5.51 

9.83 

5.22 

R 

3.56 

0.84 

3.58 

0.85 

8.74 

1.22 

8.35 

1.14 


0.30 

-0.88 

0.29 

-0.88 

5.75 

-2.10 

5.51 

-2.03 

O 

F 

-0.24 

-1.45 

-0.25 

-1.44 

6.49 

2.56 

6.21 

2.41 

p 

0.88 

-0.44 

0.85 

-0.45 

1.15 

-2.52 

1.12 

-2.43 

T 

-1.57 

-3.06 

-1.56 

-3.03 

4.52 

0.46 

4.38 

0.42 

A 

2.73 

1.24 

2.67 

1.20 

2.17 

-1.98 

2.19 

-1.91 

L 

-0.44 

-2.60 

-0.45 

-2.58 

2.01 

0.20 

1.98 

0.16 


3.47 

1.10 

3.40 

1,07 

-0.31 

-2.04 

-0.23 

-1.98 


0.30 

_^.33 

0.26 


1.48 

-0.12 

1.52 

-0.11 

A 

3.76 

0.96 

3.69 

0.94 

3.86 

1.40 

3.80 

1.37 

A 

N 

2.84 

-1.00 

2.78 

-1 .00 

0.19 

-3.04 

0.28 

-2.93 

11 

T 

3.31 

-0.38 

3.27 

-0.38 

5.15 

1.22 

4.99 

1.13 

E 

2.40 

-0.43 

2.34 

-0.46 

-0.38 

-4.51 

-0.29 

-4.34 

R 

2.12 

-0.74 

2.10 

0.73 

9.15 

4.74 

8.83 

4.50 

I 

n 

4.09 

0.34 

4.05 

0.22 

6.56 

-1.20 

6.37 

-1.19 

U 

P 

2.34 

-1.50 

2.33 

-1.49 

7.41 

0.08 

7.16 

0.03 

l\ 

2.29 

0.13 

2.27 

0.12 

4.97 
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Figure I. Reconstruction of cardiac cross section fron multi- 

planar x-ray views and by use of algebraic algorithm. 
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Figure 2. Schematic of the surgically isolated, metabol i ca I I y 
supported canine left ventricle. 
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Figure 5. Typical shape and dimension changes of a cross section during diastole, 
and finite element partitioning pattern. 
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Some node numbers identifying the four regions of extreme stresses. 
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Normalized circumferential stress distributions across the wal 
thickness at the four regions of the left ventricular cross 
sect i on . 



